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The spin Hall effect of light (SHEL) is the pho- 
tonic analogue of spin Hall effects occurring for 
charge carriers in solid-state systems. 1 3 A prime 
example of this intriguing phenomenon occurs 
when a light beam refracts at an air-glass in- 
terface. It amounts to a polarization-dependent 
displacement perpendicular to the plane of in- 
cidence. At optical wavelengths, this shift is 
about a few tens of nanometres. Recently, it was 
predicted that a light beam projected onto an 
oblique plane can undergo a significantly larger 
displacement: 4,5 This effect, named geometric 
SHEL, is a consequence of spin-orbit coupling and 
is largely independent from the physical imple- 
mentation of the projection. Here, we experi- 
mentally demonstrate this novel phenomenon by 
observing an optical beam transmitted across an 
oblique polarizer. The spatial intensity distribu- 
tion of the transmitted beam depends on the in- 
cident state of polarization and its centroid un- 
dergoes a positional displacement exceeding one 
wavelength. This novel type of spin-orbit cou- 
pling is unavoidable for light fields carrying an- 
gular momentum and practically relevant for any 
precision measurement using polarized beams. 

Already in 1943, Goos and Hanchen observed that the 
position of a light beam totally reflected from a glass- 
air interface differs from metallic reflection. 6 This is the 
most well-known example of a longitudinal beam shift 
occurring at the interface between two optical media. In 
honour of their seminal work, any such deviation from 
geometrical optics occurring in the plane of incidence, 
is referred to as a Goos-Hanchen shift. Conversely, a 
similar shift occurring in a direction perpendicular to the 
plane of incidence is known as Imbert-Fcdorov shift. 7,8 
These phenomena generally depend both on properties 
of the incident light beam and the physical interaction 
occurring at an interface. 9 Goos-Hanchen and Imbert- 
Fcdorov shifts have been observed at dielectric, 10 semi- 
conductor 11 and metal 12 interfaces. 

The Imbert-Fedorov shift 7,8,13 is an example of the 
spin-Hall effect of light (SHEL) arising at a planar in- 
terface between two dielectric media. SHEL is a conse- 
quence of the spin-orbit interaction (SOI) of light, namely 
the coupling between the spin and orbital angular mo- 
mentum in an optical field. 1,2,4,14-22 All electromagnetic 
SOI phenomena in vacuum and in locally isotropic media 
can be interpreted in terms of the geometric Berry phase 
and angular momentum dynamics. 23 For a freely propa- 
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Figure 1. Spin Hall effects of light occur when a light 
beam propagates across an oblique interface. The pres- 
ence of a surface not perpendicular to the direction of beam 
propagation breaks the symmetry of the incident beam. This 
can give rise to spin-orbit coupling, which manifests itself as 
a spin-dependent displacement of the transmitted beam. In 
this work, we specifically study light beams impinging under 
an angle 9 onto a polarizing interface. The reference frame 
{x' ,y' , z'} is aligned with the polarizer surface. 

gating paraxial beam of light, SOI effects vanish unless 
a relevant breaking of symmetry occurs. A typical ex- 
ample of such a symmetry break is the interaction of the 
beam with an oblique surface, whose normal lies under 
an angle 9 with respect to the beam axis (see Figure 1). 

Although the resulting phenomena essentially depend 
on the type of the interaction with the surface, 24 there 
are some common characteristics that reveal universality 
in SOI of light due to the geometry and the dynamical 
angular momentum aspects of the problem. Amongst 
the various observable effects resulting from the beam- 
surface interaction, the so-called geometric Hall effect of 
light is virtually independent from the properties of the 
surface 4,5,25,26 and, therefore, represents the ideal candi- 
date for studying above-mentioned universal features. It 
amounts to a shift of the centroid of the intensity distri- 
bution represented by Poynting-vector flow of the beam 
across the oblique surface of a tilted detector. Direct 
observation of this effect as originally proposed depends 
critically on the detector's response to the light field. 
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Figure 2. The geometric SHEL at a polarizing in- 
terface. (A) The horizontal dark grey plane represents 
the polarizing interface with the Cartesian reference frame 
{x',y' = y,z'} attached to it. The axis x' is taken parallel 
to the absorbing axis of the polarizing interface. 6 € [0, 7r/2) 
denotes the angle of incidence. The central wave vector fco 
of the incident beam defines the direction of the axis z of 
the frame {x,y,z} attached to the beam. It is instructive 
to study a wave vector fe in the z-y-plane, rotated by an an- 
gle 5 <C 1 with respect to fco- (B) The unit vector a _L fc, 
representing the direction of the absorbed component of the 
incident field, lies in the common plane of fc and x' (coloured 
light brown in the figure) and can be obtained, in the first- 
order approximation, from the rotation by the angle St&nO 
around fc, of the unit vector x. 



However, the question whether the response function of 
a real detector is indeed proportional to the Poynting 
vector density is subject to a long-standing debate. 27,28 

In this work, we implement an alternative scheme, 5 in 
which the occurrence of the beam shift is independent of 
the detector response. To this end, we send a circularly 
polarized beam of light across a tilted polarizing inter- 
face to demonstrate a novel kind of geometric Hall effect, 
in which the centroid of the resulting linearly polarized 
transmitted beam undergoes a spin-induced transverse 
shift up to several wavelengths. This approach is differ- 
ent from the early proposal 4 in that it is observable with 
standard optical detectors. Furthermore, it is different 
from the SHEL occurring in a beam passing through an 
air-glass interface 2 since this geometric SHEL is practi- 
cally independent of Snell's law and the Fresnel formulas 
for the interface. 

As for conventional SHEL, the physical origin of this 
geometric version resides in the SOI of light. In fact, we 
can describe this effect in terms of the geometric phase 
generated by the spin-orbit interaction as follows: When 
the beam passes through the polarizing interface, the 
plane wave components with different wave vectors ac- 
quire different geometric phases determined by distinct 
local projections yielding to effective "rotations" of the 



polarization vector around the wave vector. The interfer- 
ence of these modified plane waves produces a redistribu- 
tion of the intensity spatial profile of the beam resulting 
in a spin-dependent transverse shift of the intensity cen- 
troid. 

This effect can be better understood with the help of 
Figure 2 that illustrates the geometry of the problem. 
The incident monochromatic beam can be thought of as 
made of many plane wave components with wave vectors 
fc spreading around the central one fco = kz, which rep- 
resents the main direction of propagation of the beam, 
where |feo| = k = \k\. For a well-collimated beam, 
the angle between an arbitrary wave vector fc and the 
central one fco is, by definition, small: i$ « 1. In the 
first-order approximation with respect to <5, we consider 
k = k(z cos S — y sin S) = k(z + n y y), where n y = k y /k = 
— sin<5 = — <5. The either left- (er = +1) or right-handed 
(ct = — f) circular polarization of the incident beam is 
determined by the unit vector u a = (x + iay)/^/2 glob- 
ally defined with respect to the axis z of the beam frame 
{x,y,z}. However, from Maxwell's equations it follows 
that the divergence of the electric field of a light wave 
in vacuum is zero. This requires that the polarization 
e CT (fc) of each plane wave component of wave vector fc 
must necessarily be transverse, namely fc • e CT (fc) = 0. 
This requirement is clearly not satisfied by u„ for which 
one has fc • u a /k = ian y /\/2 ^ 0. Anyhow, the trans- 
verse nature of the light field can be easily restored by 
subtracting from u„ its longitudinal component: u a — > 
e a (k) oc u a — k (fc • u a ) /fc 2 = u a — (ia K y j\pi)z. Now 
that we have properly modeled the polarization of the 
incident field, let us see how it changes when the beam 
crosses the polarizing interface. 

A linear polarizer is an optical device that absorbs ra- 
diation polarized parallel to a given direction, say x', 
and transmits radiation polarized perpendicular to that 
direction. The electric field of each plane wave com- 
ponent of the beam sent through the polarizing inter- 
face, changes according to the projection rule e a (k) — > 
e a (k) — a(a ■ e CT (fc)) = t(t ■ e CT (fc)), where a = x — 
k(k ■ x')/k 2 is the effective absorbing axis with a = 
a/\a\ = x — yK y ta,xi9, and i = a x fc/fc is the 
effective transmitting axis. The key term to calcu- 
late is the transmitted amplitude t ■ e a (k) cx (1 — 
ian y tan Q)j\[2 = exp(— ian y teoT.9)/^/2, where an irrele- 
vant Ky-independent overall phase factor has been omit- 
ted. 

Therefore, as a result of the transmission, the am- 
plitude of each plane wave component is reduced by 
a factor 1/V2 and multiplied by the geometric phase 
term exp(— icr^j, tan 6). Here, the "tan^" behavior of the 
phase, characteristic of the geometric SHEL, 4 is in strik- 
ing contrast to the typical "cot 9" angular dependence 
of the conventional SHEL as, e. g., the Imbert-Fedorov 
shift. 29 However, the spin-orbit interaction term cx/ty, ex- 
pressing the coupling between the spin a and the trans- 
verse momentum n y of the beam, is characteristic of both 
phenomena, thus revealing their common physical origin. 



3 



The appearance of the geometric phase term 
exp(— iaKy tan 6) leads to a beam shift along the axis 
y. This can be seen by writing the incident circularly 
polarized paraxial beam as \& m (y) = u^^iy), where 
only the dependence on the relevant transverse coor- 
dinate y has been explicitly displayed, and ij)™{y) = 
J V , CT n ( K j/) cxp(ikn y y)dK y . After transmission across the 
polarizing interface, the Fourier transform of 
ip™(y) changes to ex P( — i&Ky tan0)/-\/2 and the 

output field can be written as ip° ut (y) = (1/V2) J ip^i^y) 
exp[ikn y (y — a tan6/k)]dn y = (1 / \^2)ip]^ (y — cr tan#/fc). 
This last expression clearly shows that the scalar ampli- 
tude of the output field is equal, apart from the (l/\/2) 
factor, to the amplitude of the input field trans versally 
shifted by crtan(?/fc. In summary, for circular input po- 
larization the position of the centroid of the transmitted 
beam is given by 



state preparation 

P^OC Piezo 



shift measurement 



(y) 



Jy\r, ut (y)\ 2 dy 

J |^out( y )|2 dy 



■ tan ( 
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where k = ■ (Please refer to the Supplementary Infor- 
mation for a detailed calculation of this result.) 

In the experiment, we use a Corning Polarcor polarizer, 
made of two layers of elongated and oriented silver nano- 
particles embedded in a 25 mm x 25 mm x 0.5 mm glass 
substrate. Interaction with these particles effectively po- 
larizes the transmitted beam. In order to avoid parasitic 
effects from the glass surfaces (ng = 1.517), we have sub- 
merged the polarizer in a tank with index matching liquid 
(Cargille laser liquid 5610, = 1.521). Without this liq- 
uid, the effective tilting angle inside the glass polarizer 
would be limited by Sncll's law to arcsin(l/riG) ps 41°. 

The polarizer used in this work can be described by its 
absorbing axis - the orientation of the nano-particles - 
and two empirical parameters describing transmittance 
and extinction ratio as a function of the tilting angle 9. 
Please refer to the SI for an experimental characterization 
of the device and a phenomenological model. 

In our setup (Figure 3), a fundamental Gaussian light 
beam (A = 795 nm) is prepared with the state of polariza- 
tion alternating between left- (<r = +1) and right-hand 
(a = —1) circular. By simultaneously observing the po- 
sition of the light beam transmitted across a polarizing 
interface and the incident state of polarization, we can 
determine the helicity-dependent beam shift A = 2(y) ■ 

It is possible to observe different beam shifts in two 
configurations. First, we measure the displacement At = 
^(y)\E T \ 2 °^ ^ e total transmitted energy density distri- 
bution \Et\ 2 when switching the incident state of polar- 
ization from er = +ltoer = — 1 (Figure 4(a)). Then, 
we employ an additional polarization analyzer in front of 
the detector and observe the shift A y = 2{y)^ E ^ of the 

energy density \E y \ 2 = \Et • y\ 2 of the vertically polar- 
ized field component solely (Figure 4(b)). These variants 
of the experiment coincide for polarizers with perfect ex- 
tinction ratios but can differ significantly for real-world 




Stokes r# 
measurement 



Figure 3. Experimental Setup: Simultaneous measure- 
ment of the incident state of polarization and the 
position of a light beam transmitted across a tilted 
polarizer. State preparation: The relative phase between 
horizontally and vertically polarized components is modu- 
lated using a polarizing beam splitter (PBS), a piezo mirror 
(M), quarter and half wave plates (QWP, HWP). The beam 
is spatially filtered using a single-mode fibre (SMF). Stokes 
measurement: The transmitted port of a beam splitter (BS) 
is used to monitor the state of polarization. We acquire the 
Stokes parameters S3, which distinguishes left- and right-hand 
circular polarization. Shift measurement: The beam is prop- 
agated across our sample, a tank containing a glass polarizer 
and an index-matching liquid, and its position is observed 
using a quadrant detector. An optional PBS (A) can be em- 
ployed as an analyzer in front of the detector. The photo 
currents I t , and I b are amplified and digitally sam- 

pled. 



polarizers with minor deficiencies. 

The beam shift observed in the latter case in- 
creases proportionally to the tangent of the tilting 
angle, exceeding one wavelength. This characteristic 
"tan 0" -dependence, especially close to grazing incidence 
6 90°, is unique to the geometric spin Hall effect of 
light. 

To the best of our knowledge, this is the first direct 
measurement of this intriguing phenomenon. The geo- 
metric spin Hall effect of light should not be confused 
with the conventional SHEL or Imbcrt-Fcdorov shift. 
The latter occurs at a physical interface and, while such 
interfaces are present in our experimental setup, they can 
only give rise to beam shifts significantly smaller than the 
observed effect. In Figure 5, we compare our results to 
the Imbert-Fedorov shift, which could occur at the polar- 
izer substrate, for a set of realistic parameters. 18,29 This 
illustrates that the beam shifts measured in this work 
constitute a novel spin Hall effect of Light, virtually in- 
dependent from surface effects. 

In conclusion, we have presented a novel beam shift 
occurring at an oblique polarizing interface. To test this 
geometric spin Hall effect of light experimentally, a circu- 
larly polarized laser beam was propagated across a suit- 
able polarizer. We have demonstrated that the centre 
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Figure 5. Conventional spin Hall effect of light com- 
pared to the geometric SHEL. A light beam transmit- 
ted across an interface between two media can undergo a 
transverse displacement known as the Imbert-Fedorov effect 
or conventional SHEL. Here, we plot this displacement for a 
left-hand circularly polarized beam (a = +1) for two different 
cases and compare this with the geometric SHEL studied in 
this work, (a) and (b) Geometric SHEL |At and 5 Ay for 
the two configurations studied experimentally (Figure 4). (c) 
Geometric SHEL as predicted for an ideal polarizing interface, 
(d) Conventional SHEL occurring at an air-glass interface 
(ni = 1, ri2 = 1.5). (e) Conventional SHEL expected for the 
entrance face of our submerged polarizer (m = til = 1.521, 
ri2 = tig = 1.517). 



Figure 4. Polarization-dependent beam shifts occur- 
ring at a tilted polarizer. Measurement data and two 
theoretical curves are shown. The dashed blue line shows the 
theory for an ideal polarizer, while the solid lines were cal- 
culated from a realistic polarizer model as derived in the SI. 
Error bars indicate the statistical standard error of the mean, 
(a) Overall displacement At of the intensity barycentre after 
transmission across the polarizer, (b) Displacement A y of 
the vertically polarized intensity component solely. 

of mass of the transmitted light field is displaced with 
respect to position of the incident beam as predicted by 
the theory. While a Gaussian light beam itself is invari- 
ant with respect to rotation around its axis of propaga- 
tion, the geometry induced by the tilted polarizer, breaks 
this symmetry. This amounts to a spin-to-orbit coupling 
characteristic for spin Hall effects of light. 



METHODS 

The state preparation setup (Figure 3) was aligned 
such that the intensities of vertically and horizontally 
polarized components were constant and equal. The rel- 
ative phase between these components was modulated 
with a frequency of / mo d = 29 Hz. We have used a single 
mode fibre to clean the spatial profile and two half-wave 
plates to compensate for the polarization change due to 



the fibre. Those wave plates before and after the fibre 
were aligned such that the combined operation of wave 
plates and fibre is reduced to a phase shift. We colli- 
mated the light beam using an aspheric lens (New Focus 
5724-H-B) and aligned its focus such that the beam waist 
is at the position of the detector. 

In order to simultaneously measure the beam position 
(y) and the incident state of polarization, we employed 
a dielectric mirror (Layertec 103210) as a non-polarizing 
beam splitter. At an angle of incidence of 3°, transmit- 
tance, reflectivity, and phase shift coincide for s and p 
polarization within experimental accuracy. We have ver- 
ified this important characteristic of our beam splitter 
by simultaneously measuring the states of polarization 
of the transmitted and reflected beams. To this end, we 
can replace the sample in the transmitted beam path 
with another Stokes measurement setup. 

The beam reflected from the beam splitter was propa- 
gated across the sample and its position was measured us- 
ing two vertically adjacent elements of a quadrant photo- 
diode (Silicon Sensors QD50-6). From the amplified 
photo currents It, It, the vertical position of the light 
beam can be calculated as (y) = / y+y ■ The factor / 
depends on the beam diameter and was measured in-situ 
by translating the detector using a micrometre stage. 

The transmitted port of the beam splitter was used to 
monitor the Stokes parameter S3 = I f + + I f~ ■ We identi- 
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fied S3 = 0.99 ± 0.01 and S3 = -0.99 ± 0.01 with the cir- 
cular states of polarization a + and a" respectively. The 
amplified photo currents I t , h, I a + , and l a - were dig- 
itally sampled (using a National Instruments PCI-6281) 
for Is at 50kHz. 

Since the signal is periodic with the modulation fre- 
quency /mod, we can filter technical noise in a post- 
processing stage. To this end, the discrete Fourier trans- 
form is computed and only spectral components with fre- 
quencies equal to / mo d an d higher harmonics thereof are 
passed. 

For both circular states of polarization, a = +1 and 
a = —1, we calculated the mean of all corresponding 
beam positions (y) and the helicity-dependent beam shift 
5 R = (y) (a = +1) - (y) {a = -1). Both series of shift 
measurements were repeated five times. We report the 
mean value and standard deviation of the mean in Fig- 
ure 4. 



We have extensively characterized systematic errors in 
our system and found that the observed position of the 
light beam depends on the state of polarization even if 
no sample is present in the beam path. The magnitude 
of this apparent beam shift ranges from 5 to 400 nm and 
depends on the transverse alignment of the collimating 
lens and the presence of polarization optics in the beam 
path. To account for this effect, the shift measurements 
A(0) = A R (#) — A R (0°) reported here are corrected with 
respect to the raw data points A R (#). 
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SUPPLEMENTARY INFORMATION 

A. Phenomenological polarizer model and 
characterization of our real-world polarizer 

In this section, we describe a geometric polarizer 
model, analogously to the work by Fainman and 
Shamir, 30 and determine empirical parameters relevant 
for the actual polarizer used in our experiment. Here, 
the interaction of an arbitrarily oriented polarizer with 
a plane wave is discussed while we deal with real light 
beams in the subsequent section on beam shifts. 

A polarizer is an optical device that alters the state of 
polarization and intensity of a plane wave without effect- 
ing its direction of propagation k = k/k. The polarizer 
used within this work can be described by a real- valued 
unit vector P a describing its absorbing axis. Projecting 
P a onto the transverse plane of the electric field yields to 
an effective absorbing axis 



a(k) 



P a - (P a ■ k) k 



1 - {P a ■ «) 2 



(SI.l) 



Thus, the electric field transmitted across an idealized 
polarizer is 



E T = E T - a[a - Et) = t [t ■ E 



(SI.2) 
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Supplementary Figure 1. Transmittance of light beams 
across the polarizer as a function of the tilting angle 
9 and the incident state of polarization. The polarizer 
is aligned such that at normal incidence (8 = 0°) vertical (V) 
polarization is transmitted and horizontal (H) polarization is 
blocked. We observe how the transmittance changes when the 
polarizer is rotated around the vertical axes and compare the 
experimental data (circles) to our phenomenological model 
(SI.4) (solid lines). 



B. Calculation of beam shifts occurring at an 
oblique polarizer according to our phenomenological 
model 



where Ej = Ej(k) is the amplitude of the plane wave 
exp(i k ■ r), and 



t(k) = a(k) x k 



\Pn X k\ 



(SI.3) 



is the effective transmitting axis. 

However, real-world polarizers have finite extinction 
ratios and an experimental characterization shows that 
the effectiveness of our polarizer decreases when tilted 
(Supplementary Figure 1). Thus, we have found it con- 
venient to phcnomcnologically describe the transmitted 
light field as 



E T = T t t[t- E 



a [a- Ei 



(SI.4) 



Here, we employ two empirical parameters, T a {0) and 
Tt (ff) , depending on the angle 9 between the propagation 
direction k and the unit vector perpendicular to the sur- 
face of the polarizer. We have found the following set 
of parameters to be in good agreement with both the 
observed transmission and beam shifts: 



T t {9) =1-1.2 exp(-12 cosfl) 
T a {6) = 0.51 exp(-1.3 cos (2) 



(SI.5) 
(SI.6) 



A more detailed study of tilted polarizers will be pub- 
lished elsewhere. 



In this section, we adapt the theory for the geometric 
spin Hall effect of light originally calculated for a different 
type of polarizer 5 to our phenomenological model. To 
this end, we express the polarizer's absorbing axis 

P a = x' = cos^ x + sin 6z (SI.7) 

in the global reference frame {x,y,z}, aligned with the 
direction z of beam propagation. 

The incident beam Ej{r) is circularly polarized and 
well-collimated, i.e. it has a low divergence 9q, with a 
Gaussian transverse intensity profile. This is expanded 
in a plane wave basis with amplitudes Ej(k) such that 
the /^-dependent projection (SI.4) can be applied. 

As a consequence of our polarizer model, the electric 
field Ex{r) after transmission across such a polarizer, is a 
superposition of two orthogonally polarized field compo- 
nents, E x (r) and E y (r). Thus, the electric field energy 
density distribution (at the detection plane z = 0) 



\E T (x,y)\ 2 = \E T -x\ 2 + \E T -y\ 



\E T -z\ 2 



(SI.8) 



= |£x| 2 



can be decomposed analogously. 

Geometric SHEL manifests itself as a transverse dis- 
placement (y) of the transmitted light beam's barycen- 
tre. The total electric energy density \Ex\ 2 and both 
of its non- vanishing components undergo such a shift. 
Since this spatial displacement is independent from the 



7 



z coordinate, we restrict the discussion to the detection 
plane at z = 0. It is convenient to write the total energy 
density's baryccntre 



w 



<y>, 



w y (y) 



(SI.9) 



as a weighted sum of the relative shifts occurring for the 
horizontally and vertically polarized components respec- 
tively. Here, 

JJ yu(x,y)dxdy 



{V)u 



II u(x,y)dxdy 



(SI.10) 



denotes the centre of mass along the y direction calcu- 
lated with respect to a scalar distribution u(x, y). The 
integration spans the whole detection plane. 
The weights 



H\E\*{x,y)dxdy 



=>v(*) = 



<D{6l) and 



(SI.11) 



(SI.12) 



introduced in (SI.9) depend on the empirical parameters 
Tt and r a (SI. 5). Finally, we can calculate the relative 
shifts 



(v) 



\E.\ 2 



JJy\E t \ 2 dxdy 
II 



tan# 
2tt 



fm + o(9 2 ), 



where the factors 



/«(*) 



fy(0) 



1 



n(0) 

Ta(0) 
Tq(0) 

rt{0) 



< and 



> 



(SI.13) 

(SI.14a) 
(SI.14b) 



depend critically on the performance of the polarizer. 

Note that, since the transmission coefficients are real 
and positive and 1 > r t > r a (Supplementary Figure 1), 
these relative shifts have opposite signs. Consequently, 
the displacement of the total energy density 



(y)\E T \ 2 



A 



tang 

O- — (W X f X + Wy fy) 



2ir 



(SI.15) 



<i 



is smaller than the relative shift (y)< E , 2 01 the verti- 
cally polarized component solely. For our realistic polar- 
izer model, both shifts are smaller than expected for the 
ideal case. This matches very well with our experimental 
observation. 

For a perfect polarizer with r a = and r< = 1, 
equation (SI.15) reduces to equation (1) since w x f x = 

^■^1 — — s- for r a — ^ 0. This is exactly the same ex- 
pression that was originally found for the polarizer model 
by Fainman and Shamir. 5 



C. Calculation of the geometric spin Hall effect of 
light at a polarizing interface 

In this part, we derive in a detailed manner equation 
(1) for a paraxial fundamental Gaussian beam incident 
at an angle 9 upon the surface of an absorbing polarizer 
whose absorbing axis is directed along x' . The Cartesian 
coordinate systems attached to the polarizing interface 
and to the beam are still defined as in Fig. 2A in the 
main text. However, now the axis x' is taken parallel to 
the absorbing axis of the polarizer, in order to reproduce 
the actual experimental conditions. 

Consider the fundamental solution of the scalar parax- 
ial wave equation 31 that we indicate with tp[r): 



1/2 



1 + iz/L 



exp 



1 



y 



iz/L 
(SI.16) 



where the Rayleigh range L of the beam can be expressed 
in terms of the beam waist wq as L = tog/2, with 0q = 
2/ (too) *C 1 denoting the angular spread of the beam. 
In the first-order approximation with respect to 9q, the 
electric vector field of the incident beam can be written 

oo-32 



* m (r) 



u a ip(r) + z-Ua ■ V^(r) 
. , 6»o (x + iay) 



' y/2w (1 + iz/L) 



ip(r) 



(SI.17) 



where u a = (x + iay)/y2 and an overall (irrelevant) 
multiplicative term has been omitted. In Ref. 33 it was 
demonstrated that the field transmitted by an arbitrar- 
ily oriented polarizer can be written as a perturbativc 
expansion of the form 

* out (r) = Goo* in (r-) 



Gio^^(r)+G i^^(r) 



Ox 



dy 



0{6l), 



where the 3x3 matrices G nm arc defined as 

k n+m. d n + m (nh) 



n 



\m\ dkV:dk, 



0, Aln, 



=0 



For our absorbing polarizer, the dyadic form 
nn =| n v n x ni 



IVy ll X -"y 



ri'x fly Tlx rx z 
n y n z 

2 



(Sl.lf 



(SI.19) 



(SI.20) 



n: 



is defined in terms of the effective-transmission unit vec- 
tor n = a x k = n x x + n y y + n z z, where 



a = 



x! — k(k ■ x') 
y/\-\k- x'\ 2 



(SI.21) 



8 



and k = k/k. A straightforward calculation furnishes 



for the transmitted field: 



Goo = 



Gn 




GlQ — 



tan 6 
tan 6* -1 
0-10 




(SI.22a) 



(SI.22b) 



(SI.22c) 



Substitution of equations (SI. 22) and (SI. 17) into equa- 
tion (SI. 18) yields to the following first-order expression 



* out (r) 



fir) 

V2 



xV 

yia ^1 

V 



y a tan 6 
w 1 + iz/L 

y a tan ( 



wo 1 + iz/L 



ZU 



wo 1 + iz/L 



O{0l). (SI.23) 



The electric field energy density, the physical quantity 
actually measured by a standard optical detector, is pro- 
portional to |\& out (r)| 2 which can be calculated from Eq. 
(SI.21) as: 



|* out (r)|^ = Wr)| a ( \ 



wq 1 + iz/L 



(SI.24) 



Finally, the sought shift is calculated as the first moment 
of the electric field energy density distribution, namely 



(y) 



J J y\^(r)\ 2 dxdy 



|¥ out (r)| 2 da:cto 



(7 

— tan ( 

k 



(SI.25) 



where the integration is evaluated over all the xy plane 
at z = 0. Equation (SI.25) reproduces the result given 
by equation (1) in the main text. 



